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ABSTRACT 

We report on a search with the IceCube detector for high-energy muon neu- 
trinos from GRB080319B, one of the brightest gamma-ray bursts (GRBs) ever 
observed. The firebaU model predicts that a mean of 0.1 events should be de- 
tected by IceCube for a bulk Lorentz boost of the jet of 300. In both the direct 
on-time window of 66 s and an extended window of about 300 s around the GRB, 
no excess was found above background. The 90% CL upper limit on the num- 
ber of track- like events from the GRB is 2.7, corresponding to a muon neutrino 
fluence limit of 9.5 x 10~^ergcm~^ in the energy range between 120 TeV and 
2.2 PeV, which contains 90% of the expected events. 

Subject headings: gamma rays: bursts - methods: data analysis - neutrinos - 
telescopes 



1. Introduction 

Long-duration gamma-ray bursts (GRBs) are thought to originate from the collapse of 
a massive star into a black hole thereby releasing a huge amount of energy in 7-rays into 
the surrounding medium. Assuming an isotropic emission of these 7-rays, the measured 
fluxes yield an isotropic equivalent energy of C(10^^-10^^ erg). However, the actual released 
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energy can be significantly lower if the 7-rays are only emitted within a small cone (jet) 
as suggested by the observation of signatures for jet breaks in some X-ray spectra. Apart 
from being amongst the most violent events in the universe, GRBs also belong to the few 
plausible source candidates for ultra-high-energy cosmic rays. Though our knowledge about 
GRBs has greatly increased in recent years, their exact nature, the way in which particles are 
accelerated, and the composition and generation of the jets formed from material accreted 
onto the black hole are still not fully understood. The observation of high-energy neutrinos 
from GRBs would be a smoking gun evidence for the acceleration of hadrons in the jets and 
hence for the connection between GRBs and extra-Galactic cosmic rays. 

In the fireball model (Mcszaros & Rces 1993), neutrinos of energy 0(10^^ eV) are pro- 
duced in the decay of charged pions generated in the interaction of accelerated protons 
of energy O{10^^ eV) with keV-MeV photons via the A"*" resonance (Waxman & Bahcall 
1997). Both synchrotron and inverse Compton^ emission from accelerated electrons have 
been proposed as the mechanism for the production of these photons which form the 7- 
ray signal measured by satellites. The particle acceleration is thought to occur in internal 
shocks (Narayan et al. 1992; Rees & Meszaros 1994; Sari & Piran 1997) yielding spectra 
for protons and electrons as typically expected in Fermi acceleration (Waxman & Bahcall 
1997). The energy in protons (normalization of the proton spectrum) is usually quoted in 
relation to the energy in electrons which is linked to the energy in 7-ray photons through 
the synchrotron and inverse Compton energy-loss mechanisms. In the pion decay neutrinos 
are produced with the flavor ratios (1/^:1/^:1/^) — (1:2:0) ^, changing to (1:1:1) at the Earth 
due to oscillations^. Their fluence is proportional to the fluence in 7-rays. 

The first calculations of the expected neutrino fiux from GRBs (Waxman & Bahcall 
1997; Alvarez-Muniz & Halzen 1999) used average GRB parameters and the GRB rate mea- 
sured by BATSE to determine an all-sky neutrino fiux from the GRB population. This 
so-called Waxman-Bahcall GRB fiux or similar GRB fiuxes have been probed with the 
AMANDA-II neutrino telescope (Achterberg et al. 2007, 2008) with negative results. These 
fiuxes will be detectable by next-generation neutrino telescopes like IccCubc with an instru- 
mented volume of > 1 km^. However, the average fiux for a single burst derived in these 



"'^Production of 7-rays through inverse Compton emission for bursts with low optical luminosity (majority 
of bursts) is actually disfavored according to Piran et al. (2009) as it leads to a very high energy 7-ray 
component in the TeV range which would carry much more energy than the observed prompt 7-ray emission 
resulting in an "energy crisis" with most current progenitor models. 

^Here and throughout the rest of the paper, u denotes both neutrinos and antineutrinos. 

^Kashti & Waxman (2005) showed that above a certain energy (typically 0(10 PeV)) the ratio changes 
to (1:1.8:1.8) due to cooling energy losses of the muons producing the neutrinos. 
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models is very small even for km^ detectors. Nevertheless, as the expected neutrino flux can 
actually vary by orders of magnitude between GRBs due to fluctuations in the burst param- 
eters ( Alvarez- Muniz et al. 2000; Becker et al. 2006), the detection of extremely bright GRBs 
in neutrinos does seem possible, albeit requiring at least a km'^-size detector, as, for example, 
was demonstrated in the analysis of GRB 030329 with the AMANDA detector (Stamatikos 
et al. 2005). 

On 2008 March 19, at 06:12:49 UT GRB080319B (Racusin et al. 2008a) was detected by 
the Swift (Burrows et al. 2005) and Konus-Wind (NASA 1994) satellites at R.A. = 217? 9 and 
decl. = 36? 3. It was the optically brightest GRB ever observed and with a peak magnitude 
of 5.3 even visible to the naked eye for a short period of time, despite the fact that the 
corresponding redshift was about 0.9. It is also one of best measured GRBs with optical 
wide-field observations covering the whole duration of the explosion (Cwiok et al. 2008) and 
with many triggered follow-up observations spanning the electromagnetic energy spectrum 
from radio to 7-rays (Racusin et al. (2008a) and references therein). 

2. Neutrino spectrum calculation 

We calculate the expected prompt neutrino spectrum for GRB 080319B in the internal 
shock scenario of the fireball model following Guetta et al. (2004) which is based on Waxman 
& Bahcall (1997). This model allows for the easy incorporation of many measured parameters 
of GRB080319B (in Waxman & Bahcall (1997) average GRB parameters are used) which is 
crucial when investigating a GRB that deviates largely from the average GRB. Other models 
like Murase & Nagataki (2006)^ or Razzaque et al. (2003) do not provide this possibility 
without obtaining the actual simulation code and are therefore not considered here. 

For reference, we list all formulae used in our calculations in Appendix A. One of the 
major inputs to the model is the keV-MeV 7-ray spectrum recorded by the satellites. In 
contrast to Guetta et al. (2004) where a broken power law was used, we parameterize the 
7-ray spectrum with a Band function (Band et al. (1993); Equation (A2)). The function 
parameters obtained from a fit to the time integrated Konus-Wind spectrum are taken 
from Racusin et al. (2008a) (suppl. information) and are hsted in Table 1. The table also 
contains additional parameters with their values required by the model. Not all of them are 
measured or even well known. For the jet parameters (fraction of jet energy in electrons), 



^This model is actually similar to Waxman & Bahcall (1997) but uses Monte Carlo simulation to calculate 
the photomeson production in interactions and takes the synchrotron losses of mesons and protons into 
account. 
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Ev (GeV) 

Fig. 1. — Calculated muon neutrino spectrum for different bulk Lorentz boost factors Fjct of 
the GRB jet: Tjct = 300 (solid line), Tjet = 500 (dash-dotted hue), and Tjct = 1400 (dotted 
line). For comparison, the average Waxman-Bahcall GRB fluence for a single GRB is also 
shown (fine-dotted line). 

eb (fraction of jet energy in magnetic field) and /e (ratio between energy in electrons and 
protons) typical values of 0.1 are used (Becker 2008). The variabihty of the 7-ray flux iyar, 
which is used as a measure for the time between the emission of two consecutive shells, was 
analyzed in Margutti et al. (2008). They find an initial time scale of 0.1 s which increased 
to 0.7 s in the course of the emission. However, their analysis was limited to the energy 
range between 15 and 150 keV, where for high energies (100-150 keV) the dominant time 
scale was significantly shorter (0.05s). By contrast, muons reconstructed in the IceCube 
detector are mainly produced by neutrinos near the first break energy (Figures 1 and 5(b)) 
which originate from proton interactions with 7-rays around 500 keV (break energy in 7-ray 
spectrum). Due to the large gap between 150 and 500 keV and the lack of an extrapolation 
method we use a typical value of tyar = 0.01 s (Becker 2008) in our calculations. 

The neutrino spectrum is parameterized as a Band function with a broken power law at 
high energies (Equation (A5)), where the latter describes the steepening of the spectrum due 
to synchrotron losses of pions and muons. The energy in neutrinos, i.e., the normalization of 
the spectrum, is proportional to and the energy in 7-ray photons (Equation (AlO)). The 

Fjct factor, which enters the model equations to the second and fourth power (Equations (A6), 
(A7), and (A9)), has a large impact on the normalization of the neutrino fiuence. With 
increasing Fjct, the shell collisions occur at larger distances from the black hole where the 
photon field and hence the target density for the pion production is smaller. The value of 
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Fjet can be estimated from the fact that the source has to be transparent for 7-rays near the 
maximum 7-ray energy produced (Guctta et al. 2004). Using the parameters given in Table 1 
and a maximum 7-ray energy of 100 MeV yields Fjct ~ 300. In Racusin et al. (2008a), the 
authors argue that the exceptional brightness of the optical flash in GRB 0803 19B implies 
that the self-absorption frequency cannot be far above the optical band and they determine 
Fjet to lie between 300 and 1400. Pandey et al. (2009) obtain Fjet ~ 300 from an extrapolation 
of the late-time evolution of the afterglow. Using a synchrotron self-Compton model, Kumar 
& Panaitescu (2008) find a Fjet factor of ~ 500. For the calculation of the neutrino spectrum, 
we adopt the optimistic case with Fjet = 300 which is displayed in Figure 1 as the solid line. 
An increase of the Fjet factor to 500 (1400) decreases the neutrino fiuence by about a factor 
10 (10^) and shifts it to higher energies (Figure 1 dashed and dotted lines). Muon cooling 
(Kashti & Waxman 2005) affects the neutrino spectrum only at energies above ~ 20 PeV 
(Fjet = 300) and is therefore negligible for our analysis. 



3. Analysis of IceCube data 

IceCube (Achterberg et al. 2006), the successor of the AMANDA experiment and the 
first next-generation neutrino telescope, is currently being installed in the deep ice at the 
geographic South Pole. Its final configuration will instrument a volume of about 1 km^ of 
clear ice in depths between 1450 m and 2450 m. Neutrinos are reconstructed by detecting the 
Cherenkov light from charged secondary particles, which are produced in interactions of the 
neutrinos with the nuclei in the ice or bedrock. The optical sensors consist of photomultipliers 
housed in pressure- resistant glass spheres (digital optical modules (DOMs); Abbasi et al. 
(2009)) which are mounted on vertical strings. Each string carries 60 DOMs with the final 
detector containing 80 such strings. Physics data taking with IceCube started in 2006 with 
nine strings installed. At the beginning of 2007, the detector was enlarged to 22 strings. 
Since 2009 April, IceCube has been running with 59 strings. The completion of the detector 
construction is planned for the year 2011. 



3.1. Data sets, event reconstruction and selection 

The data acquisition (DAQ) system of IceCube (Abbasi et al. 2009) is based on local 
coincidences between neighboring DOMs in a string within 1 /xs for which the photon signal 
passes a threshold of 0.4 photo electrons. All data from DOMs belonging to a local coinci- 
dence are read out and the digitized waveforms are sent to a computer farm at the surface. 
In order to pass the trigger, a minimum number of eight DOMs in local coincidences within 



Table 1. Fireball Model Parameters Used in the Calculation of the Neutrino Spectrum for 

GRB 080319B 



Parameter Value Reference 







1.3 X 10^4 erg 


Racusin et al. (2008a) 


Burst duration 


66 s 


Racusin ct al. {2008b) 


Tjet 




300, 500, 1400 


Racusin et al. (2008a); sec also main text 


7 spec. 


— fluence, (20keV-7MeV) 


6.23 X 10"*crgcm-2 


Racusin et al. (2008a) 


7 spec. 


— breaJj energy, e-y 


651 keV 


Racusin et al. (2008a) (suppl. information) 


7 spec. 


— 1st index, aj 


0.83 


Racusin et al. (2008a) (suppl. information) 


7 spec. 


— 2nd index, (3^ 


3.5 


Racusin et al. (2008a) (suppl. information) 


z 




0.94 


Vreeswijk et al. (2008) 


T 




0.2 


Becker (2008) 


Ee 




0.1 


Becker (2008) 


es 




0.1 


Becker (2008) 


fe 




0.1 


Becker (2008) 


tvar 




0.01s 


Becker (2008) 



Fraction of proton energy going into pion in a single p7 interaction. 



Table 2. Number of Expected Signal and Background Events at Different Cnt Levels 







Signal 


Background (Off-Time Data) 




Cut Level 


No. Events 


Efficiency'' (%) 


No. Events'' Efficiency^ (%) 


Comment 


Trigger 


0.24 


100 


1.2 X 10-1 100 


See Section 3.1 


Quality 


0.14 


58 


1.4 X 10-3 1.2 


See Equation (1) 


Final 


0.10 


41 


1.7 X 10-5 0.014 


See Equation (2) 



'■Relative to trigger level. 

^In a cone with radius 5° centered on GRB080319B position within 66 s. 



-10- 




5 -0.25 — I I I I I I I 'i I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I 



-250 -200 -150 -100 -50 50 100 150 200 250 
Seconds since BAT Trigger Time 

Fig. 2.— 7-ray emission from GRB080319B as measured by Swift-BAT (Stamatikos 2008) 
with Tq — 06:12:49 UT. The dashed vertical hnes mark the time range covered by IceCube 
data. 



a time window of 5000 ns is required. If this condition is fulfilled, the waveforms are com- 
bined to an event and the number and arrival times of the Cherenkov photons are extracted. 
For each event, an initial track is reconstructed using the line- fit algorithm (Ahrens et al. 
2004). This is a simple but fast track reconstruction based on the measured hit times in the 
DOMs. We do not consider events where only a shower is produced (e.g., in interactions 
of electron neutrinos or neutral current neutrino interactions) as only the track-like light 
pattern of muons allows for a good angular resolution. Hence, this search focuses on muon 
neutrinos from GRB080319B. Cuts on the reconstructed zenith angle (> 70°) and number 
of hit DOMs (> 10) reject downgoing atmospheric muons and reduce the event rate down to 
the 117 Hz. This allows to apply more advanced reconstruction algorithms to the remaining 
neutrino candidates. A more precise estimate of the direction of a neutrino candidate is 
obtained by fitting a muon-track hypothesis to the pattern of the recorded Cherenkov hght 
in the detector using a log-hkelihood reconstruction method (Ahrens et al. 2004). A fit of 
a paraboloid to the region around the minimum of the log-likelihood function yields an es- 
timate of the uncertainty on the reconstructed direction. The absolute time of an event is 
determined with a GPS clock with a precision of better than a millisecond, which is more 
than sufficient for this analysis. 

At the time of GRB080319B, IceCube was running in maintenance mode with 9 out 
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Fig. 3. — (a) Difference between the event rate in a 5s bin, R, and the average rate R 
(calculated from all events in the 300 s time window shown) as a function of time at trigger 
level. The solid and dashed vertical lines mark the satellite trigger time of the GRB and the 
start and stop times of the measured 7-ray emission, respectively, (b) Histogram of {R — R) 
divided by the statistical errors. The hne is a fit of a Gaussian distribution to the histogram. 



of 22 strings taking data. Apart from the reduced number of strings the DAQ system had 
a slightly different configuration than during normal operation. IceCube data are available 

in a window of about 300 s around the GRB (on-time data) as displayed in FigTirc 2. The 
detector was checked for stability during this period by plotting the rate of events passing 
the trigger in bins of time, R, relative to the average rate R of 117 Hz (Figure 3). The 
variations in the event rate arc compatible with statistical fluctuations and there were no 
indications for abnormal behavior of the detector during the period under consideration. 

In order to avoid systematic uncertainties due to inaccuracies in the simulation when 
calculating the significance of a deviation from the background-only hypothesis, the expected 
background in the on-time window is determined from the observed off-time data. However, 
the amount of data taken with the special DAQ configuration during the regular mainte- 
nance runs is limited and not sufficient for a good background estimation. Instead, we utilize 
the IceCube data set of the 2006 data taking period (131 days lifetime), when the detector 
consisted only of the same nine strings which were taking data during GRB080319B. Af- 
ter applying the following quality cuts, the 2006 data set shows good agreement with the 
GRB 080319B maintenance data set (Figure 4) 



^rec > 90° 



/9.O for 4 < A^dir < 10 

(Tdir < 10 ; 6'min > 70 ; Lred < \ 

[8.0 for A^dir < 4 



with 



• ^rec- reconstructed zenith angle; 

• (7dir: uncertainty on the reconstructed track direction (quadratic average of the minor 
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Fig. 4. — Comparison after quality cuts (Equation (1)) of the background data set (2006 
IceCube 9 strings) with a 1 hr data set taken 1 week after ORB 080319B with equivalent 
DAQ settings (2008 IceCube 9 strings). 



and major axis of the la error ellipse); 

• Lred' — logio of the likelihood value of the reconstructed track divided by the number of 
degrees of freedom (number of hit DOMs minus number of fit parameters) . In conjunction 
with the selection of upgoing tracks it has proven to be an efficient variable for separating 
upgoing atmospheric neutrinos from misreconstructed downgoing atmospheric muons. It 
exploits the fact that for a light pattern originating from a downgoing muon the incorrect 
upgoing track hypothesis yields rather low likelihood values; 

• iVdir: number of photons detected within a —15 to +75 ns time window with respect to 
the expected arrival time for unscattered photons from the muon-track hypothesis; and 

• ^min: minimum of the two zenith angles from a fit of a two-track hypothesis to the light 
pattern. 

The ^min cut rejects events where two downgoing muons from independent atmospheric 
showers pass through the detector in quick succession and produce a light pattern that fakes 
an upgoing track. The difference in the overall rates of the two samples is about 10% which 
is within the statistical error of the total number of events in the maintenance data set. 

For the final analysis, we use the method described in Section 3.3 to maximize the 
potential for a discovery: the cuts on Udir and iVdir are tightened until we reach a maximum 
in the probability to detect the fluence calculated in Section 2 (Fjct = 300) with a significance 
of at least Aa. In order to avoid biasing the results, only the off-time data are used and the 
on-time data are kept "blind" during this procedure. The optimized cuts are 

f9-0 for 8 < A^dir < 10 

^rec > 90° ; CTdir < 5° ; Omin > 70° ; L,ed < < - cl _ . 2 

(8.0 for Ndir < 8 
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Fig. 5. — (a) Effective area at final cut level for muon neutrinos from the direction of 
GRB080319B as a function of the neutrino energy, (b) Expected mean number of Monte 
Carlo signal events after final cuts as a function of neutrino energy. The two vertical lines 
mark the central interval containing 90% of the events. 

After these cuts the data sample with an event rate of about 5 x 10~^ Hz is still dominated 
by misreconstructed downgoing muons (a Monte Carlo simulation of atmospheric neutrinos 
yields an event rate of 2 x 10~^Hz). For a search cone with 5° radius centered on the GRB 
position and a time window of 66 s, a mean number of background events of 1.7 x 10~^ is 
expected. A summary of the event rates at different cut levels is given in Table 2. 



3.2. Monte Carlo 

Signal muon neutrinos from GRB080319B are generated with a Monte Carlo and 
weighted according to the fluence calculated in Section 2. The Monte Carlo contains a 
detailed simulation of the propagation of the muon neutrino through the Earth and the 

ice using the ANIS generator (Gazizov & Kowalski 2005). After the interaction, the muon 
is traced through rock and ice taking into account continuous and stochastic energy losses 
(Chirkin & Rhode 2004). The photon signal in the DOMs is determined from a detailed 
simulation (Lundberg et al. 2007) of the propagation of Cherenkov light through the ice 
which includes the modeling of the changes in absorption and scattering length with depth 
(Ackermann et al. 2006). This is followed by a simulation of the DOM electronics and the 
trigger. The DOM signals are then processed in the same way as the real data. Background 
events do not have to be simulated as these are taken from off-time data. 

After final cuts (Equation (2)) and for neutrinos from the direction of the GRB (weighted 
according to the calculated GRB spectrum), 90% of all reconstructed Monte Carlo tracks 
are contained within 20° of the true direction. The median angular resolution is 5.6°. This 
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resolution is worse than the one usually quoted for IceCube in its nine-string configuration for 
neutrino-induced muons. The reason for this is that the geometry of the detector was such 
that the reconstruction lever arm was at its shortest for the direction of GRB080319B. The 
corresponding muon neutrino effective function of energy is displayed in Figure 5(a). 

The expected mean number of events from the GRB (Fjet = 300) after final cuts is 0.1, with 
90% of the events contained in the energy range from 120 TeV to 2.2 PeV (Figure 5(b)). A 
summary of the event rates at different cut levels is given in Table 2. 

In order to compare the signal (neutrino) Monte Carlo to data, a high-purity (atmo- 
spheric) muon neutrino sample from the off-time data is selected by requiring zenith > 100°, 
(Jdir < 2? 5, and A^^ir > 8 (these tight cuts are only used for this comparison and not in 
the following analysis of the data). Figure 6 shows good agreement between data and the 
Monte Carlo weighted to the atmospheric muon neutrino fiux of Barr et al. (2004) (Bartol 
spectrum) . 

3.3. Unbinned likelihood analysis 

The final analysis is based on data sets at final cut level (Equation (2)). The data are 
analyzed using an unbinned log-likelihood method similar to the one described in Braun 
et al. (2008). In contrast to binned methods where the event is rejected if it lies outside the 
cut region (binary selection), unbinned likelihood methods do not throw away events but 
use probability density functions (PDFs) to evaluate the probability of an event to belong 
to signal or background. Therefore, they are more powerful than binned methods. 

In the case of searches for neutrinos from GRBs detected by sateUites, the direction 
and time of reconstructed tracks are the most crucial information to separate signal from 
background. Hence, both the signal, S{xi), and background, B{xi), PDFs are each the 
product of a time PDF and a directional PDF, where Xi denotes both the directional and 
time variables. The directional signal PDF is a two-dimensional Gaussian distribution with 
the two widths being the major and minor axis of the la error ellipse of the paraboloid fit 
described in the previous section. The time PDF is fiat over the 7-ray emission time and 
falls off on both sides with a Gaussian distribution with cr = 25 s. The Gaussian accounts 
for possible small shifts in the neutrino emission time with respect to that of the 7-rays and 
prevents discontinuities in the likelihood function. The sensitivity of the analysis depends 
only weakly on the exact choice of a, e.g., the quoted upper limit changes by less than 2% 
when doubling a. For the directional background PDF, the detector asymmetries in zenith 
and azimuth have to be taken into account. This is accomplished by evaluating the data 
in the detector coordinate system. The directional background PDF is hence derived from 
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Fig. 6. — Comparison between an atmospheric muon neutrino Monte Carlo (solid line) and 
a high-purity data set of (atmospheric) muon neutrinos (filled squares) in different variables. 
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the distribution of all background events after final cuts in the zenith-azimuth plane of the 
detector. The time distribution of background tracks during the GRB can be assumed to be 
constant resulting in a flat time PDF. 

Both PDFs are combined in an extended log-hkelihood function (Barlow 1989) 

N 

ln(£(K))) = -{ris) - {rib) + J]ln((n,) Six^) + {n^) B{xi)) , (3) 

i=l 

where the sum runs over all reconstructed tracks left after cuts with Xi representing the 
PDF parameters (absolute time of the track along with the track direction in local detector 
coordinates and its estimated uncertainty). The variable {nj,) is the expected mean number 
of background events, which is determined from the background data set. The mean number 
of signal events, (n^), is a free parameter which is varied to maximize the expression 

in order to obtain the best estimate for the number of signal events, (fis)- 

To determine whether a given data set is compatible with the background-only hypoth- 
esis, a large number of background data sets for the on-time window are generated from the 
2006 nine-string data by randomizing the track times while taking into account the down- 
time of the detector. For each of these data sets, the ln(7^) value is calculated, yielding 
the distribution shown in Figure 7(a). The probability for a data set to be compatible with 
background is given by the fraction of background data sets with an equal or larger ln(7^) 
value (p- value). The sensitivity of the IceCube detector to neutrinos from GRB080319B is 
determined by injecting different numbers of Monte Carlo signal events into the generated 
randomized background data sets and calculating the resulting p- value (see Figure 7(a)). 

GRBs are expected to generate a substantially harder neutrino energy spectrum than 
that of atmospheric neutrinos. A detector quantity closely related to the neutrino energy 
is the number of DOMs (channels) with detected photons, Nchan- This quantity is used 
to enhance the sensitivity to a possible signal. Figure 7(b) shows the A^chan distribution 
of off-time data events within a cone with 20° radius centered on the GRB position (no 
events remain in the on-time data set) together with a signal Monte Garlo which has been 
normalized to the number of data events. Combining this information with the p- value from 
the ln(7^) distribution increases the detection probability. The potential for discovering a 
muon neutrino fiuence as calculated in Section 2 (Fjet = 300) with a significance of ia or 
larger is about 6%. 
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Fig. 7. — (a) Likelihood ratio distributions for background-only (off-time) data sets (solid 
line) and for the same background data sets with one (dashed line) and two (dotted line) 
signal events injected, (b) Nchan distributions for off-time data (squares), and signal Monte 
Carlo (solid line) after final cuts. Data events are selected within a radius of 20° around the 
GRB position. The Monte Carlo has been normalized to the number of events in data. 

3.4. Results and systematic uncertainties 

After all parameters of the analysis have been fixed, the on-time data are unblinded and 
analyzed with respect to two time windows^. The shorter one (from Tq — 3.8 s to Tq + 62.2 s) 
corresponds to the immediate emission time of the 7-rays, whereas the second (from Tq — 173 s 
to Tq + 130 s) covers the whole time range with IceCube data (see Figure 2). No significant 
deviation from the background-only hypothesis was found in either of the two time windows. 
In both cases, the unbinned likehhood method yields ln(7^((n5))) — with {ris) = as 
the best estimate for the number of signal events. The resulting Neyman 90% CL upper 
limit (Neyman 1937; Amsler et al. 2008) on the number of signal events in the short time 
window is 2.7, i.e., injecting signal events according to a Poisson distribution with mean 2.7 
into the on-time data yields \n(JZ{{ns))) > in 90% of cases. From the neutrino fluence 
calculated in Section 2 for Fjct = 300 we expect to see 0.1 events. This results in a model 
rejection factor (MRF; Hill & Rawlins (2003)) of 27 which is defined as the ratio between 
the upper limit on the number of signal events and the expected number of signal events. 
The corresponding upper limit on the muon neutrino fluence is plotted in Figure 8(a) in 
the energy range containing 90% of the detected signal events. Integrating the fluence over 
this energy range yields an upper limit of 9.5 x 10~^ergcm~^. This upper limit is slightly 
conservative as we do not consider the effect of u-,- from GRBs. Tau neutrinos might manifest 



^The time windows define the flat part of the signal time PDFs. 
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Fig. 8. — (a) Ninety percent CL upper limit on the fluence from GRB080319B (solid line) 
with respect to the calculated neutrino fluence for a bulk Lorentz boost of the jet of Fjct = 300 
(dashed line), (b) Ratio between upper limit and calculated fluence spectrum (MRF) as a 
function of Fjet- 

themselves as r tracks (which can travel a substantial distance at PeV energies) or as muons 
from tau decays^. The impact of a larger bulk Lorentz boost on the MRF is displayed in 
Figure 8(b). In addition, the MRFs for Fjet = 500 and 1400 are listed in Table 3. 

The effects of systematic errors on the result, described below in detail, were investigated 
by varying simulation parameters and repeating the full analysis. The quadratic sum of all 
systematic errors is (+17%,— 4%). The main sources of systematic uncertainty are as follows. 

Ice simulation: Inaccuracies in the ice simulation can lead to a wrong estimate of the 
efficiency of the detector to neutrinos from the GRB. In order to estimate the size of this 
effect, the analysis was repeated using a modified ice simulation. In this simulation, the DOM 
efficiency was altered as a function of depth according to the differences observed between 
data and Monte Carlo in the DOM occupancy, effectively making the ice clearer for depths 
with above-average transparencies and dirtier for depths with below-average transparencies. 
This leads to an increase of the fluence upper limit of 16%, which is included as a one-sided 
systematic error. 

DOM efficiency: Uncertainties in the efficiency of the optical modules in the photon 
detection lead to an uncertainty in the number of expected events from the GRB. Varying 
the efficiency by ±10% changes the upper hmit by ±4%. 

Background rate: Even after optimized cuts, the data set is dominated by misrecon- 
structed downgoing atmospheric muons. The rate varies throughout the year due to changes 



^Electron neutrinos do not contribute to the signal as the resulting electron immediately produces an 
electromagnetic cascade. 
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in the density profile of the atmosphere at high altitude above the South Pole by about 10% 
around the mean value. As the number of events after cuts in the data set during the GRB is 
too low, the rate at the time of the GRB was determined using a 1 hr data set recorded with 
the same DAQ settings about a week later. In order to account for potential differences, the 
background data rate was varied by ±10%. This results in a shift of the upper limit of less 
than ±1%. 



4. Comparison with other results 

In Murase (2008), the author calculates neutrino spectra for a burst like GRB080319B 
with different models. The expected flucnces in the best cases are below our optimistic 
scenario. Therefore, our limits do not constraint these models. 

In Thrane et al. (2009), the Supcr-Kamiokande collaboration reports on upper limits on 
the neutrino fluence from GRB080319B. They analyzed their data in different energy ranges 
between ~100MeV and ~lTcV. Their strongest fluence limits come from the high-energy 
upmu data set and amount to 16cm~^ (22 cm^^) at 90% CL for (anti) muon neutrinos as- 
suming an E^"^ spectrum. In comparison, IceCubc reaches an upper limit of 1.1 x 10^^ cm^^ 
in the energy range E^, — 120TeV-2.2 PeV. A direct comparison of the results is difficult as 
the Super-Kamiokande paper does not contain information about the neutrino energy inte- 
gration range used in the calculation of the limits. Assuming a lower integration boundary of 
-E*!/,! = 10 GcV and an upper boundary E'^, 2 ^ E^^i yields a sensitivity for an E'"^ spectrum 
that is about 200 times worse than that of IceCube (for equal fluxes of muon and anti muon 
neutrinos) . 

Along with high-energy neutrinos which origin from the decay of charged pions high- 
energy 7-ray photons are produced in the decay of simultaneously generated neutral pions. 
In addition, high-energy photons are produced in inverse Compton scattering of synchrotron 
photons by accelerated electrons (Falcone et al. 2008). In contrast to neutrinos, the flux 
of high-energy 7-ray photons at the Earth is signiflcantly reduced due to the large optical 
depths for photon-photon pair production inside the source for not too large jet Lorentz 
factors Fjet ^ 800 (Falcone et al. 2008). In addition, high-energy photons above 100 GeV are 
absorbed on the extragalactic background light if they travel distances z > 0.5. Observa- 
tions with air-Cherenkov telescopes like H.E.S.S. (Aharonian et al. 2009; Tam et al. 2009) or 
MAGIC (Albert et al. 2007; Galante et al. 2009) arc also hampered by the fact that usually 
it takes more than 50 s (MAGIC) or 100 s (H.E.S.S.) from the observation of a GRB by a 
satellite to the start of data taking with these telescopes. Therefore, the prompt emission 
window is only partially covered or not at all. Milagro (Atkins et al. 2004; Abdo et al. 2007) 
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as an air shower array observed the visible sky continuously. However, it was mostly sensi- 
tive to energies above 1 TeV. Up to now, there has been no definitive detection of very high 
energy 7-ray emission from GRBs. Milagrito (Atkins et al. 2000, 2003) and the HEGRA 
AIROBICC array (Padilla et al. 1998) reported evidence at the 3a level for high-energy 
7-ray emission from GRB970417A {E^ > 650 GeV) and GRB920925C {E^ > 20 TeV), re- 
spectively. However, subsequent searches for high-energy 7-ray emission from GRBs did not 
find similar signals. For GRB 080319B, no high-energy 7-ray data exist (MAGIC, which has 
a threshold below 100 GeV, was not able to observe the burst as it was already dawning 
(Lorenz 2008)). Current fiux predictions are near or below the sensitivity of current instru- 
ments (Falcone et al. 2008), where the predicted fiuxes in the energy range below ~ 100 TeV 
are dominated by the leptonic emission component in most scenarios. Therefore, the upper 
limits from high-energy 7-ray observations do not constrain the neutrino fiux in our model 
or the computed upper limit. 



5. Conclusions and outlook 

We used the IceCube neutrino telescope to search for high-energy muon neutrinos from 
GRB080319B, one of the most spectacular and well measured GRBs ever observed. Based 
on the fireball model of GRBs and the measured 7-ray fiuence we calculated the expected 
neutrino spectrum for different jet bulk Lorentz boosts Fjet- After applying quality cuts 
to suppress misreconstructed atmospheric muons a mean number of 0.1 signal events is 
expected for the optimistic case of Fjet = 300 (for other Fjet see Table 3) in IceCube, which 
was running in a nine-string configuration. The data were analyzed with an unbinned log- 
likelihood method utilizing the directional and temporal distance of reconstructed tracks 
to the GRB. The sensitivity to a potential GRB signal was enhanced by evaluating energy 
information. No deviation from the background-only hypothesis was found either in a small 
time window covering the immediate 7-ray emission time or an extended window of 300 s. 
This results in a 90% CL upper limit on the muon neutrino fiuence from GRB080319B 
within the short time window of 9.5 x 10~^ergcm~^ in the energy range between 120 TcV 
and 2.2 PeV which contains 90% of the expected signal events. The corresponding ratio 
between the upper limit and the calculated GRB spectrum (MRF) is 27. Its stabihty with 
respect to systematic uncertainties in the analysis is estimated to be (-1-17%,— 4%). The 
upper limit does not allow us to impose constraints on GRB parameters within the fireball 
model. 

In its final configuration with 80 strings the expected number of detected events in 
IceCube from a bright GRB like GRB080319B is C(l), rendering the individual analysis 
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of these rare GRBs highly interesting also in the future. Using the large number of GRBs 

observed by the Swift and Fermi (formerly GLAST) satellites, the growing IceCube detector 
will also soon be able to probe the Waxman-Bahcall or similar GRB fluxes and in the case 
of a non-detection set stringent limits. 
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Table 3. Neutrino Spectrum Parameters According to the Fireball Model for 
GRB 080319B for Different Bulk Lorentz Factors of the Jet Together with the Expected 
Mean Number of Events in the Detector and the Model Rejection Factor Obtained from 

the Analysis 



Parameter \ bulk Lorentz factor 


Tjct = 300 


Tjet = 500 


Tjet = 1400 


1st break energy, ei 


260 TeV 


710 TeV 


5.6 PeV 


2nd break energy, £2 


3PeV 


23PeV 


1.4 EeV 


1st index, 


0.5 


0.5 


0.5 


2nd index, Pi, 


2.17 


2.17 


2.17 


3rd index, 7^ 


4.17 


4.17 


4.17 


Fluence at 1st break energy 


7.3 X 10-1^ GeV-icm-2 


1.4 X 10-i*GeV-icm-2 


2.9 X 10-1* GeV^^cm-^ 


Expected events in IceCube 


0.10 


8.6 X 10"^ 


9.6 X 10"** 


MRF 


27 


420 


3.7 X 10^ 
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A. Equations used in the calculation of the neutrino spectrum 
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The parameters of the 7-ray spectrum Fy{Ey) are the break energy Cy and the spectrum 
indices before and after the break ay and (5y, respectively. The quantity Ty is the measured 
7-ray fiuence and z the redshift of the GRB. The parameters of the neutrino spectrum 
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Fi,{Ei,) are the two break energies, ei and 62- and the corresponding spectral indices aj,, 
(3^^ and 'y^. The expression 1 — (1 — {xp^Tr))^^^^^^ in Equation (AlO) estimates the overall 
fraction of the proton energy going into pions from the size of the shock, AR, and the mean 
free path of a proton for photomeson interactions, Ap^. Here, (xp^j^) = 0.2 is the average 
fraction of proton energy transferred to a pion in a single interaction. The expression ensures 
that the transferred energy fraction is < 1. The variables and f,y are obtained from the 
integrals of Equations (A3) and (AlO). The isotropic equivalent luminosity, is given by 
the isotropic equivalent energy released in 7-rays, E^°, divided by the burst duration. The 
calculations are insensitive to the beaming effect caused by a narrow opening angle of the 
jet (0?4 for GRB080319B according to Racusin et al. (2008a)) as all formulae contain the 
isotropic luminosity in conjunction with a Att shell geometry, i.e., effectively use luminosity 
per steradian. For example, the target photon density used to calculate Xp-y is given by 
oc L'^747ri?2. 
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